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Abstract. Deep observation (~ 113 hrs) of the Cygnus region at TeV energies using the HEGRA stereoscopic system of air 

^ ' Cerenkov telescopes has serendipitously revealed a signal positionally inside the core of the OB association Cygnus 0B2, at 

' the edge of the 95% error circle of the EGRET source 3EG J2033+4118, and ~ 0.5° north of Cyg X-3. The source centre 

C/3 ! of gravity is RA aj2ooo: 20'" 32" 07' ±9.2;;,^, ± 2.2;;^,, Dec Snaoa- +41° 30' 30" +2.0^t„t + 0.4;!^,. The source is steady, has a 

' post-trial significance of +4.6o", indication for extension with radius 5.6' at the ~ 3cr level, and has a differential power-law flux 

^ I with hard photon index of -1.9 ± 0.3stat ± 0.3sys. The integral flux above 1 TeV amounts ~3% that of the Crab. No counterpart 

• ^ . for the TeV source at other wavelengths is presently identified, and its extension would disfavour an exclusive pulsar or AGN 

' origin. If associated with Cygnus 0B2, this dense concentration of young, massive stars provides an environment conducive 

M , to multi-TeV particle acceleration and likely subsequent interaction with a nearby gas cloud. Alternatively, one could envisage 

■ - - I y-ray production via a jet-driven termination shock. 
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1. Introduction 

The current generation of ground-based imaging atmospheric 
Cerenkov telescopes offer coverage of the multi GeV to TeV 
y-ray sky at centi-Crab sensitivity and arc-minute resolution. 
Stereoscopy employed by the HEGRA CT-System at La Palma 
(Daum et al. 1997) offers highly accurate reconstruction of 



event directions at angles up to ~ 3° off-axis. Results here are 
taken from data originally devoted to Cyg X-3, and the EGRET 



source GeV J2035+4214 (Lamb & Macomb |1997[ ). The sepa- 
ration between these objects (~ 1.5°) permits a combined anal- 
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ysis given the overlap in their CT-System fields of view (FOV). 
This letter presents analysis details and observational proper- 
ties of a serendipitously discovered TeV source in these data. 
A brief discussion concerning astrophysical origin and location 
of this new source is also presented. 

2. Data Analysis and Results 

The HEGRA system of imaging atmospheric Cerenkov tele- 
scopes (lACT-System), consists of 5 identical telescopes oper- 
ating in coincidence for the stereoscopic detection of air show- 
ers induced by primary y-rays in the atmosphere. In data dedi- 
cated to Cyg X-3, alternate ~20 minute runs targeting the Cyg 
X-3 position +0.5° in declination were taken during moonless 
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nights of Aug-Sept 1999, Sept-Oct 2000 and Jun-Oct 2001. 
Likewise in data dedicated to GeV J2035+4124, ~20 minute 
runs were obtained tracking directly the GeV source during 
Jul- Aug 2001. In total, three tracking positions are present in 
combined data. After quality checks, a total of 112.9 hours 
data are available for analysis. Preferential selection of y-ray- 
like events (against the cosmic -ray background) is achieved 
by using the difference between the reconstructed and as- 
sumed event direction, 0, and the mean-scaled-width param- 
eter, w (Konopelko 1995). In searching for weak point-like and 
marginally extended sources, so-called tight cuts are consid- 
ered optimal given the angular resolution of the CT-System 
(< 0.1°): 6 < 0.12° and w < 1.1, where we use algorithm 
'3' as described by Hofmann et al. ( 1999 ) for the event direc- 
tion reconstruction. The number of images per event, Mtei, used 
for calculating and w was also a priori chosen at Wtei ^ 3. 
Monte Carlo simulations (Konopelko et al. 1999 ) and tests on 
real sources have shown that «tei = 2 events contribute little to 
the overall sensitivity. 



2.1. Source Search and Background Estimates 

In searching for new TeV sources, skymaps of event direction 
excesses over the RA & Dec plane are generated after having 
estimated the background over the FOV. A new empirically- 
based template background model has been developed with 
the goal of simple generation of skymaps. The template back- 
ground comprises events normally rejected according to the w 
criterion. We define the number of events in the y-ray regime 
s from w < 1.1, and for the template background b from 
1.3 < w < 1.5. A necessary correction applied to the tem- 
plate background accounts for differences in radial profile be- 
tween the two w regimes. A normalisation a, to derive excess 
events s - ab al some position in the FOV, accounts for dif- 
ferences in the total number of events in the two w regimes. 
A ful l description of the template model appears in Rowell 
2002. Fig. |l] presents the resulting excess skymap. The tem- 
plate model was used in discovering the TeV source which 
is evident ~ 0.5° north of Cyg X-3. An event-by-event cen- 
tre of gravity (COG) calculation (Table |l]a), weighting events 
with +1 from the s and ab regimes respectively is performed. 
The COG accuracy is li mited b y a systematic pointing error of 
~ 25" (Piihlhofer et al. 
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1997). A pre-trial significance at the 



COG position of +5.9o- is obtained, summing events within 
= 0.12° (Table |l]b). Statistical trial factors arise from the 
initial discovery skymap (different to that in Fig. [l|) in which 
event directions are independently summed in 1 100 bins of size 
0.1° X 0.1°. Assuming 1100 trials are accrued in locating the 
COG, the post-trial probability = 1 .0 - (1 - P)i'°° for P the 
pre-trial probability (one-sided P - 1.9 x 10"^, or -1-5. 9cr), is 
then calculated as P, = 2.1 x 10"^. This gives a post-trial sig- 
nificance of -H4.6cr. 1 100 is actually a slightly conservative trial 
estimate since oversampling of the y-ray point spread function 
(PSF) by a factor ~ 1.5 occurs in the discovery skymap. To 
verify results using the template model, we make use of a con- 
ventional type of background model employing background re- 
gions displaced from the on-source region spatially in the FOV 
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Fig.l. Skymap (1.5° x 1.5° view at 0.05° x 0.05° binning) of 
excess events s - ab, using the template background model. At 
each bin, the excess is estimated from events within a radius 
- 0.12°. Included are 95% eiTor ellipses of variou s EGRET 
sources, the core of Cygnus OB2 (Knodlseder 2000 ), the TeV 
COG (star), its 2cr error circle, and the location of Cyg X-3. 
ASCA GIS contours (2-10 keV) are overlayed. 



but derived from the same m> < 1.1 regime. Background events 
are taken from ring-segments with matching trigger character- 
istics to that of the source region. A normalisation a according 
to the solid angle ratio between background and on-source re- 
gions is then applied. Results using this so-called ring model 
(Table. |l]b) are consistent with those from the template model. 

2.2. Observational Properties of the TeV Source 

Splitting data firstly according to their three tracking posi- 
tions reveals commensurate source contributions (Table |l|c). 
The source is also found to develop linearly with the cumu- 
lative number of background events. Such tests suggest consis- 
tency with a steady source during the three years of data col- 
lection. We have also verified that after cuts a constant back- 
ground acceptance throughout the dataset is observed and that 
the event excess in vv-space appears consistent with that of a 
true y-ray population. To determine source size, we fit a ra- 
dial Gaussian convolved with the point spread function (deter- 
mined from Crab data) to the excess events as a function of 
0-^, using a subset of events with the best angular resolution 
(wtei = 5) for which errors are minimised (Fig. The intrin- 
sic size of the TeV source is estimated at ctsjc - 5.6' + 1.7'. 
Correlations between the fit parameters suggest that the signif- 
icance for a non-zero source size is at the ~ 3. Oct level rather 
than the 3.5cr level indicated above. A breakdown of the ex- 
cess with «tei also shows that the ntei=5 exclusive subset con- 
tributes strongly to the excess (Table |l]d). Such behaviour is 
suggestive of a generally hard spectral index given that higher 
trigger multiphcities are favoured by higher energy events. For 
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(b) Tight cuts: < 0.12% w < 1.1, n,ei > 3 
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(c) Tigilt cuts on tracliing subsets 



Back. 


t 


r\ s b 


a 


5 






Cyg X-3 6 - 0.5° — 






Template 


39.5 


0.69 148 647 


0.172 


+3.0 


Ring 




148 1994 


0.057 


+3.3 






- Cyg X-3 (5+0.5° — 






Template 


45.4 


1.00 276 1214 


0.170 


+4.2 


Ring 




276 1266 


0.168 


+3.8 






~~ GeV J2035 ~ 






Template 


28.0 


0.68 99 472 


0.156 


+2.6 


Ring 




99 1193 


0.057 


+3.4 



f. Observation time (hrs) 
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(d) Tight cuts on iix^ subsets 
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(e) Spectral Cutst: 6 < 0.224°, w < 1.1, n,ei > 3 
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' : Aharonian et al. |1999j summarise other spectral cuts. 

Table 1. Summary of numerical results for the TeV source, un- 
der two background models. Here, s and b are the resulting 
event numbers for the y-ray-like and background w regimes 
respectively, and s - abis the derived excess using a normali- 
sation a. S denotes the excess significance using Eq. 17 of Li 
& Ma |l983i See section |for definitions of 6 and w. 



the energy spectrum calculation and selection cuts, we follow 



the method of Aharonian et al. (1999 1 using effective collec- 
tion areas appropriate for on-axis and -1° off'-axis sources as 
per the exposure efficiency for the TeV source in these data. 
A tight cut w <1.1, as opposed to the less-restrictive w < 1.2 
is also used. For energies below ~0.8 TeV the effective collect- 
ing area decreases markedly sharper at positions beyond ~ 1 .0° 
off-axis compared to positions nearer on-axis. Limiting our fit 
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Fig. 2. Distribution of 6^ for 9 calculated from the COG (solid 
histogram & filled dots), against a background derived from the 
template model (dashed hist). The convolved radial Gaussian fit 
F ^ P3+P2 exp{-d^/(Pl^ + o-2j)) is indicated by the soHd fine 
with PI - o-,,-c =0.094° + 0.026° the intrinsic source size. The 



PSF width crpt=0.070° 
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is estimated from Crab data. 
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Fig. 3. Differential energy fluxes of the TeV source and other 
results. 'H- A' is t he AIROBICC 90% confidence level upper 
limit (Prahl |1999| ) at the TeV COG converted to diflFerential 
form at 20.8 TeV assuming a spectral photon index of -2.0. 
We interpret the 3EG J2033+4118 flux as an upper limit. The 
ASCA GIS 99% upper limit assumes a photon index of 2, and 
Nh - 10"^ cm"^. Data fits (solid lines) and their Icr statistical 
eiTors in photon index (dashed lines) are shown. TeV data are 
wen fit (jcl = 0.23/3) by a power law (Eq. |I]). 

therefore to energies >0.8 TeV reduces systematic errors. A so- 
called loose cut 6 < 0.224° using the ring background model 
is used in deriving the energy bin-by-bin excess since a loose 
cut in 6 improves the y-ray selection efficiency according to the 
moderately-extended nature of the source. Results are shown in 
Table [l^ and Fig. H, with the spectrum being well fit by a pure 
power law with generally hard photon index. Systematic eiTors 
are estimated from changes in bin centres and uncertainties in 
Monte Carlo-derived collection areas: 



phcm^^s^'TeV"^ 



(1) 



1-13 



dN/dE = N (E/neVy 

N = 4.7 (+2.1 Stat + 1.3sys) X 10- 

r = -1.9(+0.3stat±0.3sys) 

The integral flux at F(E > 1 TeV) = 4.5 (±1.3stat) x 10"'^ ph 
cm"^ s ' amounts to 2.6% that of the Crab. 

3. Discussion & Conclusion 

The OB association Cygnus 0B2 is unique for its com- 
pact nature and extreme number of member OB and O stars 
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(eg. Knodlseder pOOOl ), and both theoretical and observational 
grounds for non-thermal particle acceleration have long-been 
discusse d (eg . Montmerle |1979|, Ca sse & Paul |1980t Volk & 
Forman 1982, White & Chen 1992). The TeV source is posi- 
tioned inside the core of Cygnus OB2 as defined by Knodlseder 
(2000). Assuming the TeV source is as distant as Cygnus OB2 



(1.7 kpc), a luminosity ~ 10 erg s" above 1 TeV is implied, 
well within the kinetic energy (KE) budget of Cygnus OB2 
estimated recently by Lozinskaya et al. (2002) at a fewxlO-''' 
erg s 



' , and also within the KE budget of a number of notable 

Manchanda et 



member stars (eg. Massey & Thompson 1991 
al 



1996 Benaglia et al. 2001). So far no counterparts at other 



wavelengths are identified. No massive or luminous Cygnus 
OB2 star of note discussed recently (eg. Massey & Thompson 
1991L Romero et al. |1999[ Herrero et al. |200li Benaglia et al. 



2001) is positioned within the Icr TeV error circle. No cat- 
alogued X-ray source from the ROSAT all-sky and pointed 
survey lies within the 2cr TeV error circle. Our analysis of 
archival ASCA GIS data yields a 99% upper limit (2-10 keV) 
of 2.0x10"'^ erg cm"^ s ' (Fig. ||). Such results may imply 
that the energy source for particle acceleration is not co-located 
with the TeV source, arising instead from the winds of the 
young/massive stars of Cygnus 0B2, either individually or col- 
lectively, or from an alternative source. The former scenario 
would generally favour accelerated hadrons interacting with a 
local, dense gas cloud, giving rise to 7r°-decay TeV emission. 
The likely hard TeV spectrum can be explained by (energy- 
dependent) diffusion effects, accelerator age, and accelerator 
to target distance (see eg. Aharonian 2001). There is how- 



ever at present no strong indication from CO and HII surveys 
(Leung & Thaddeus 1992, Dame et al. 2001, Comeron & Torra 
2001) for any co-located dense gas cloud, although see Marti 
et al. (200C) who discuss a nearby HII cloud in the context 



that follows immediately below. A suggested alternative sce- 
nario involves a jet-driven termination shock at which acceler- 
ated electrons produce synchrotron a nd TeV inverse-Compton 
(IC) emission (Aharonian & Atoyan 1998 ). Such a jet could 
emanate from a nearby microquasar, possi bly a class of high 
energy y-ray source (see eg. Paredes et al. 2000 ). In fact two 
nearby sources, 3EG J2033+41 18 and also the E GRET source 
possibly associated with Cyg X-3 (Mori et al. 1997) could 



be GeV indicators of such a microquasar Rem arkably, Cyg 
X-3 appears to have a bi-lobal jet (Martf et al. |2000t |2001| ) 
well-aligned with the TeV source, the latter which would be 
~ 70pc from Cyg X-3 in absolute terms if it is at the same 
distance (>8.5 kpc). Future X-ray observations will be a cru- 
cial constraint on the IC emission in this context. We interpret 
the flux from 3EG J2033+4118 presently as a MeV/GeV up- 
per limit at the TeV COG (Fig. |^, and note that the directly 
extrapolated energy flux from the TeV source lies about two 
orders of magnitude below the 3EG J2033-I-41 18 flux at over- 
lapping energies. Finally we note that earlier claims for a TeV 
source (Neshpor et al. 1995, at a ~1 Crab flux level) and flaring 



episodes coincident with a Cyg X-3 radio flare at energies >40 
TeV (Merck 1993 Krawczynski 1995| ) positionally consistent 
with our TeV COG have been reported. These results are how- 
ever in conflict with our estimates of the flux level and steady 
nature of the TeV source assuming they all have the same ori- 



gin. Further observations with the HEGRA CT-System aimed 
at confirmation and improving our spectral and source mor- 
phology studies are now underway. 
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